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ABSTRACT 


"  Progressive  gravity  wave*  at  the  interface  between  two  unbounded  fluids 
are  considered.  The  flow  in  each  fluid  is  taken  to  be  potential  flow.  The 
problem  is  converted  into  a  set  of  int agio- differential  equations,  reduced  to 
a  set  of  algebraic  equations  by  discretisation,  and  solved  by  Newton's  aethod 
together  with  parameter  variation.  Meiron  and  Saffman's  calculations  showing 
the  existence  of  overhanging  waves  are  confirmed.  However,  the  present 
calculations  do  not  support  Saffaan  and  Then's  conjecture  that  the  waves  are 
geometrically  limited  (i.e.  that  solutions  exist  until  the  interface 
intersects  itself).  It  is  proposed  that  along  a  solution  branch  starting  with 
sinusoidal  waves  of  small  amplitude,  one  reaches  solutions  with  vertical 
streamlines  and  then  overhanging  waves.  Continuing  on  this  branch  one  returns 
to  nonoverhanging  waves  and  thence  back  toward  a  wave  with  vertical  stream¬ 
lines.  It  is  suggested  that  this  succession  of  patterns  and  accompanying 
oscillation  in  wave  characteristics  is  repeated  indefinitely.  Graphs  of  the 
results  are  included.. 
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SIGNIFICANCE  AND  EXPLANATION 


Y  ' 


Gravity  waves  propagating  on  the  interface  between  two  fluids  of 
different  densities  are  observed  in  laboratory  experiments  and  in  oceano¬ 
graphic  studies*  Internal  waves  propagating  along  theraoclines  in  the  ocean 
and  along  the  interface  between  salt  water  and  fresh  water  in  estuaries  are 
but  two  examples.  This  report  is  part  of  an  ongoing  theoretical  and  numerical 
study  of  this  phenomenon. 

Here  we  do  a  numerical  study  of  the  interfacial  periodic  gravity  waves 
progressing  along  the  interface  between  two  unbounded  fluids  of  different 
densities*  Assuming  two-dimensional  potential  flow  in  each  fluid,  the 
resulting  problem  is  reduced  to  a  set  of  algebraic  equations  by  discretisation 
and  solved  by  Newton's  method  together  with  variation  of  wave  parameters. 
Calculations  of  Vanden-Broeck  and  others  are  extended  to  show  that  progression 
along  a  one  parameter  family  of  waves  toward  a  limiting  configuration  is 
accompanied  by  oscillatory  behavior  in  the  wave  characteristics. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
lies  with  MIC,  and  ndt  with  the  authors  of  this  report. 
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OH  THZ  LIMITING  CCHFIGUJiATIOH  Or 
INTIRFACIAL  GRAVITY  MAVIS 


f.I.l.  Turner  and  J.-M  Vandan-Broack 

1.  Introduction 

«t  consider  two-dls*n«lonal  progressive  gravity  mvh  propagating  at  tha  interface 
between  two  unboondad  fluid*  of  diffarant  danaitias  p1  and  P2  where  P1  ia  tha 
danaity  of  tha  lower  fluid.  For  pj  ■  0  tha  problem  raducaa  to  that  of  tha  elaaaical 
Stokes'  wav a*,  stokaa  propoaad  that  thaaa  waves  raaeh  thair  Uniting  configuration  whan 
tha  valoeity  in  tha  horiaontal  diraetion  equal*  tha  phaaa  apaad  of  tha  wav*.  In  a 

frame  of  rafaranca  Moving  with  tha  wave,  tha  valoeity  at  tha  craat  ia  than  aqual  to  aaro 
and  tha  fra*  aurfaea  aakaa  a  120*  angla  with  itaalf.  Stokaa'  conjactura  haa  baan 
confirmed  in  tha  work  of  Mick,  rraankal,  and  Toland1.  For  p2  >  0,  Stokaa*  limiting 
configuration  is  no  longer  poaaibla  because  such  a  stagnation  point  in  tha  lower  fluid 
would  result  in  an  infinite  velocity  ia  tha  upper  fluid  at  that  point. 

Conjactura*  about  tha  limiting  configuration  of  intarfaclal  waves  war*  made  by 
Holysr3  and  Saffman  and  Yuan3.  Holyar2  performed  axtanalva  computations  by  using 
Pad*  approximanta  and  found  waves  for  which  tha  fra*  surface  profile  la  vertical  at  sosm 
point.  In  a  similar  problem  with  tha  fluids  confined  between  horiaontal  walla  it  ia  shown 
analytically  by  Mick  and  Turner4  that  if  tha  solitary  waves  do  not  become  infinitely 
broad  in  the  horiaontal  direction  than  vertical  streamline*  must  appear.  Based  on  her 
calculations  Holyar2  conjectured  that  a  wav*  with  a  vertical  profile  at  son*  point 
constituted  a  limiting  configuration.  Sine*  the  horiaontal  velocity  at  such  a  point  would 
•goal  the  phase  apaad,  Holyar's2  criterion  raducaa  to  Stokes'  criterion  when  p2  ■  0. 

and  Yuan3  pointed  out  that  thara  are  no  dynamical  or  kinanatical  reason*  to  reject 
intarfaclal  waves  for  which  tha  horiaontal  velocity  exceeds  tha  phase  velocity  at  some 


Sponsored  by  tha  united  States  Army  under  Contract  Ho.  BAAG29-80-C-0041  and  the  National 
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point.  The  existence  of  such  waves  was  aubaaquantly  demonstrated  numerically  by  Matron 
and  Saffnan5,  for  thaaa  wavaa  portion*  of  the  haaviar  fluid  lia  abova  lightar  fluid  and 
for  thia  raaaon  Matron  and  Saffaan5  callad  than  “overhanging"  wavaa  (cf  figure  3c). 

The  computations  of  Matron  and  Saffaan5  have  thus  shorn  that  the  appearance  of  a 
vertical  tangent  on  the  interface  does  not  correspond  to  a  limiting  configuration. 

Saffswn  and  Yuan3  proposed  that  the  overhanging  wavaa  continue  to  exist  until  the 
interface  intersects  itself.  However,  the  calculations  of  Matron  and  Saffaan5  ware 
stopped  well  before  any  such  Halting  configuration  and  thus  do  not  confirm  the  conjecture 
of  Saffaan  and  Yuen3. 

In  the  present  paper  we  compute  inter facial  wavaa  by  using  Vanden-Broeck*  a®  schema 
(hereafter  referred  to  as  V-B).  The  problem  is  first  formulated  as  an  integro- 
differential  aet  of  equations.  These  equations  are  than  reduced  to  a  set  of  algebraic 
equations  through  discretisation  and  the  resulting  aquations  solved  by  Mewton'a  method. 

Our  results  confirm  the  existence  of  overhanging  wavaa  as  computed  by  Nelron  and 
Saffauui5.  Moreover,  our  scheme  enabled  us  to  extend  their  calculations.  Continuing  along 
the  branch  of  solutions  we  compute,  we  find  that  there  is  a  return  from  the  overhanging 
configuration  to  one  in  which  the  fluid  interface  is  a  single-valued  function  of  the 
horizontal  variable.  Ne  believe  that  further  progress  along  the  branch  will  produce 
steepening  of  the  streamline  and  developstent  of  a  second  overhanging  wave.  Accompanying 
this  oscillation  in  shape  is  an  oscillation  in  speed  and  amplitude  as  dearly  shown  in 
Fig.  2.  He  conjecture  that  an  alternation  between  nonoverhanging  and  overhanging  waves 
continues  indefinitely  as  one  proceeds  along  the  solution  branch.  Such  oscillatory 
behavior  has  been  exhibited  in  the  analytical  work  of  longuet-Higglne  and  Fox7  on  the 


surface  wave  problem, 


2.  Formal* tlon 


let  us  oonaldsr  mvu  *t  the  interface  between  two  fluids  of  infinite  axtent  having 
danaitiaa  snd  p2,  where  is  tha  dansity  of  tha  lover  fluid*  Assume  a  wave  of 

wavelength  1  propagates  with  phase  velocity  e  undar  tha  inf luanca  of  gravity  g.  Tha 
variables  are  aada  diaansionlass  by  taking  X  as  tha  unit  length  and  c  as  tha  unit 
velocity.  Wa  choose  a  Cartesian  fraae  of  reference  in  which  tha  flow  is  steady.  The  x- 
axis  is  parallel  to  tha  velocity  at  an  infinita  distance  frost  tha  interface  and  the  y-axis 
is  directed  vertically  upward  so  that  gravity  acts  in  the  negative  y  direction.  It  is 
assessed  that  tbs  interface  is  symmatric  with  rsspect  to  tha  y  axis. 

Lst  and  +2  **•  potential  functions  and  let  and  g2  be  stream  functions  in 
the  lower  and  upper  fluids  respectively.  Nithout  loss  of  generality  we  choose 

■  ♦j  ■  0  at  a  crest  and  -  $2  ■  0  at  the  interface.  The  wave  speed  parameter, 
the  density  parameter,  and  the  steepness  are  defined  by  the  relations 


s  -  y(0)  -  y(y) 


(1) 

(2) 

(3) 


By  choosing  an  origin  for  y  so 
lower  fluids  to  the  upper  is  (P2 
at  the  interface  is  express'!  by 


that  the  difference  in  the  total  Bernoulli  head  from  the 

2 

c 

~  P1 )  — j,  the  condition  that  the  pressure  be  continuous 


~  pqjj)  +  <!  "  P>jr  *  4^  o  -  p) 


(4) 


Here  and  q2  are  the  magnitudes  of  the  dimensionless  velocities  where  “  0  and 
♦2  “  0.  Mote  that  the  choice  of  the  origin  for  y  is  such  that  in  an  undisturbed  fluid 
the  Interface  is  at  y  *  0. 

Mow  consider  x  and  y  as  functions  of  and  in  ths  lower  fluid  and  as 

functions  of  +2  and  t2  in  the  upper  fluid.  This  provides  two  different  paraarntric 
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representations  for  the  interface;  namely,  x(^ ,0-),  y(^ ,0-)  and  x(*2 ,0+),  y(+2,0+) 
where  0-  and  0+  indicate  that  tha  second  variable  has  b«an  allowed  to  approch  aaro 
from  below  and  above,  respectively.  Ha  usa  the  notation  x^  (^ ) ,  y^^)  and 
x2<«2),  y2<+2) ,  raapaetivaly,  for  thaaa  boundary  values,  in  tha  present  context  tha 
Hilbert  transform  provides  tha  integral  relations 

<**,  v2  <jy«  .  ,  ... 


-  ’  - 1  v  i~*  «♦;  -  ",'K 

■  1  *  Z4  — ?  [«*  *1*1  -  ♦„>  -  ~t  .<♦'  <■  »  >]*' 


o  d* 


(cf  V-H°). 


In  general,  ♦1  *  f2  for  two  points  in  contact  on  tha  interface.  To  relate  the 
values  of  the  potentials  we  define  the  fuqption 


♦,  -  *<♦,> 


by  the  "contact  relations" 


x2(g2)  -  x^g^));  y2<*2)  -  y,(g(*2))  (8) 

and  normalise  g  by  requiring  g(0)  -  0.  In  the  paper  V-B  the  basic  equations  were 
written  using  ♦1  as  the  independent  variable  and  a  change  of  variables  was  used  to  have 
mesh  points  equally  spaced  in  +2,  thereby  concentrating  points  where  the  velocity  in  the 
upper  fluid  changes  rapidly.  Hare  the  same  end  is  accomplished  by  using  +2  as  the  basic 
independent  variable.  By  differentiating  the  relations  in  (8)  one  can  rewrite  equation 


47  •  P]^]2  +  U  -  PW  ’  ft  °  ’  P> 


and  aquation  (5)  aa 

dx  V2  dy 

d*J  ^df^"1  "  1  ■  f  ^4  (cot  w  [si<*2>  ■  g<*2^  ■  oot  *  +  ,(*2,la*2  (10) 


Th*  problH  ctn  now  be  formulated  as  follow*  t  fix  two  of  th*  throe  parameter* 

P,  V,  and  at  than  find  function*  x2(f2),  Ya (♦j > »  end  g($2)  defined  for 
♦2  *  1.0,  l^J,  and  a  value  for  the  remaining  parameter  ao  that  (6),  (9),  and  (10)  are 
simultaneously  satisfied.  We  solve  the  problem  numerically.  The  relations  (6),  (9)  and 
(10)  are  replaced  by  algebraic  equations  to  be  satisfied  at  discrete,  equally  spaced 
values  of  th*  variable  end  th*  resulting  algebraic  equations  are  solved  using 

Newton's  method.  Tor  more  details  regarding  the  numerical  procedure  see  V-B.  In  that 
paper  p  and  s  were  given  and  V  was  found  as  part  of  th*  solution.  To  obtain  th* 
results  presented  her*  we  used  the  scheme  from  V-B  as  well  as  a  scheme  with  th*  additional 
option  of  fixing  p  and  u  and  finding  a  as  part  of  the  solution. 


3.  numerical  results 

The  scheae  outlined  in  the  previous  section  was  used  in  V-B  to  coapute  interfacial 
waves  for  p  «  0.1  with  H  <  25,  where  N  denotes  the  nuaber  of  mesh  points  between 
♦2  *  0  and  As  can  be  seen  from  table  1,  the  results  in  V-B6  agree  with  those 

2  i 

of  Holyer  .  Saffaan  and  Yuen  repeated  the  calculations  by  using  a  numerical  scheme  based 
on  series  truncation  and  table  1  shows  agreement  to  five  decimal  places  with  the  results 
of  Holyer2  and  V-B6  providing  a  valuable  check  on  all  three  procedures. 


s 

£ 

Holyer 

Saffaan  &  Yuen 

Vanden-Broeck 

(H  -  25) 

N  -  30 

N  -  60 

eP 

2* 

.880628 

.880626 

.880624 

.880625 

.880626 

.8 

21 

.930496 

.930484 

.930484 

.930484 

.930484 

Table  1 .  Selected  values  of  the  speed  parameter  P  as  a  function  of  steepness 
s  for  p  «  0.1  from  Saffaan  and  Yuen3,  Vanden-Broeck6,  and  present 
calculations.  Saffaan  and  Yuen3  use  H  ”  2ws  and  [  ( 1  +  p)|i/(1  -  p)]  ^2 
for  steepness  and  speed,  respectively. 

In  V-B  solutions  were  obtained  starting  from  saall  amplitude  waves  and  continuing  up 
to  a  wave  with  a  profile  which  is  almost  vertical  at  sene  distance  frea  the  crest  (cf. 

V-B,  figure  1).  Overhanging  waves  were  not  obtained  in  V-B  due  to  the  use  of  relatively 
few  mesh  points.  The  computations  froa  V-B6  are  continued  in  the  present  work,  wherein 
more  mesh  points  are  used  (N  -  30,  45,  60)  and  continuation  along  a  branch  of  solutions  is 
achieved  by  switching  between  u  and  s  for  parameter isation.  The  convergence  with 
increasing  values  of  N  can  be  seen  froa  the  selected  valves  in  table  1. 

Our  results  with  H  -  60  agree  with  those  of  Meiron  and  Saffaan6  and  confirm  the 
existence  of  overhanging  waves.  Moreover,  our  scheae  enabled  us  to  carry  forward  Meiron 
and  Saffaan* s5  calculations.  Fig.  1  shows  the  solution  pairs  (s,u)  for  p  taking  the 
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values  0.1,  0.S,  and  0.7.  The  broken  line  (in  rig.  2)  corresponds  to  t'.ie  results  obtained 
by  Matron  and  Saffman5  for  p  “  0.1  while  the  solid  curve  shows  our  extension  of  their 
results.  Fig.  2,  for  p  -  0.1,  shows  in  soma  detail  the  oscillation  referred  to  in  the 
introduction  and  will  be  discussed  more  fully  in  the  following  paragraphs. 

Apart  from  the  shift  in  the  solution  set  with  changing  p,  which  can  be  seen  in  Fig. 
1,  the  characteristics  of  the  waves  at  various  points  along  each  solution  branch  were 
substantially  unaltered  by  changes  in  p.  To  envision  the  variation  in  wave 
characteristics  for  a  fixed  p  it  is  useful  to  envision  moving  on  a  solution  branch  using 
arclength  as  a  parameter.  we  describe  the  behavior  for  the  case  p  ■  0.1  which  appears 
typical  and  for  which  we  have  the  most  extensive  data. 

Starting  with  an  infinitesimal  steepness  a  and  a  value  of  V  ”  ~  “  .818 

progression  along  the  branch  initially  produces  an  increase  in  both  s  and  p  (cf  figure 
1,  p  ■  0.1)  a  typical  wave  profile  along  this  section  is  shown  in  figure  3a,  the  letter 
"a"  corresponding  to  the  point  "a"  in  figure  2.  The  increase  in  s  and  v  on  this  first 
section  is  accompanied  by  an  increase  in  the  auxiaram  slope  of  the  Interface  up  to  the 
point  "b"  at  which  the  interface  has  a  vertical  tangent.  The  corresponding  wave  is  shown 
in  figure  3b.  Continuing  past  "b*  one  reaches  a  maximum  for  s  and  shortly  thereafter  a 
maximum  for  |i.  Fast  point  "b"  and  up  to  the  point  "e“  the  wave  profile  is  no  longer  a 
graph  over  x,  but  is  "overhanging"  as  pictured  in  figure  3c  corresponding  to  point  "c" 
where  there  appears  to  be  the  greateet  overhang.  At  point  "d"  there  is  a  slight  overhang 
and  at  point  "e"  it  appears  to  have  disappeared  so  that  the  profile  at  "e"  resembles  that 
at  "b".  In  progressing  past  point  "e"  one  sees  further  evidence  of  oscillation  in  s  and 
U«  The  computed  profiles  past  "e"  had  no  overhang  and  resembled  the  wave  in  figure  3a. 
Fast  the  last  point  shown  in  figure  2  the  computations  became  unreliable  and  were 
discontinued. 
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Conclusion 

The  results  of  Longue t- Higgins  and  Fox7  show  that  the  speed  of  surface  waves 
(corresponding  to  p  <■  0)  oscillates  infinitely  often  as  the  wave  of  greatest  height  is 
approached.  Here  ws  have  given  numerical  evidence  that  a  similar  oscillation  occurs  for 
.'.nterfacial  waves  with  density  ratio  p  >  0.  whereas  the  steepness  did  not  oscillate  in 
ths  calculations  of  Longuet-Higgins  and  Fox7,  we  find  evidence  of  oscillation  in  both  the 
speed  and  steepness  parmeters  in  the  case  that  p  >  0.  The  accompanying  profiles  of  waves 
which  go  from  nonoverhanging  to  overhanging  and  thence  back  to  the  nonoverhanging 
configuration  lead  us  to  conjecture  that  the  limiting  behavior  is  characterised  by 
infinitely  many  oscillations  in  all  of  the  characteristics  of  the  waves. 
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Progressive  gravity  waves  at  the  interface  between  two  unbounded  fluids  are 
considered.  The  flow  in  each  fluid  is  taken  to  be  potential  flow.  The  problem 
is  converted  into  a  set  of  in teg ro- differential  equations,  reduced  to  a  set  of 
algebraic  equations  by  discretisation,  and  solved  by  Newton's  method  together 
with  parameter  variation.  Meiron  and  Saffman's  calculations  showing  the 
existence  of  overhanging  waves  are  confirmed.  However,  the  present  calculations 
do  not  support  Saffman  and  Yuen's  conjecture  that  the  waves  jure  geometrically 


ABSTRACT  (cent.) 


20.  limited  (i.e.  that  solutions  exist  until  tha  interface  intersects 
itself) .  It  is  proposed  that  along  a  solution  branch  starting  with 
sinusoidal  waves  of  small  amplitude ,  one  reaches  solutions  with 
vertical  stre anilines  and  then  overhanging  waves.  Continuing  an  this 
branch  one  returns  to  nonoverhanging  waves  and  thence  back  toward  a 
wave  with  vertical  streamlines.  It  is  suggested  that  this  succession 
of  patterns  and  accompanying  oscillation  in  wave  characteristics  is 
repeated  indefinitely.  Graphs  of  the  results  are  included. 
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